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When an electron moves in a smoothly varying non-collinear magnetic struc- 
ture, its spin-orientation adapts constantly, thereby inducing forces that act on 
both the magnetic structure and the electron. These forces may be described by 
electric and magnetic fields of an emergent electrodynamics [1—4]. The topolog- 
ically quantized winding number of so-called skyrmions, i.e., certain magnetic 
whirls, discovered recently in chiral magnets [5-7] are theoretically predicted to 
induce exactly one quantum of emergent magnetic flux per skyrmion. A mov- 
ing skyrmion is therefore expected to induce an emergent electric field following 
Faraday's law of induction, which inherits this topological quantization [8]. Here 
we report Hall effect measurements, which establish quantitatively the predicted 
emergent electrodynamics. This allows to obtain quantitative evidence of the 
depinning of skyrmions from impurities at ultra-low current densities of only 
10 6 Am~ 2 and their subsequent motion. The combination of exceptionally small 
current densities and simple transport measurements offers fundamental insights 
into the connection between emergent and real electrodynamics of skyrmions in 
chiral magnets, and promises to be important for applications in the long-term. 

Skyrmion lattice phases (SLPs) in chiral magnets such as MnSi and other B20 transition 
metal compounds are a new form of magnetic order, composed of topologically protected 
vortex lines (the skyrmions) with a non-zero winding number that are stabilized parallel 
to a small applied magnetic field. Skyrmion lattices in magnetic materials were discovered 
only recently by means of small angle neutron scattering (SANS) [5, 9]. The winding num- 
ber of the spin structure was thereby first inferred from topological contributions to the 
Hall effect [6] and Lorentz force microscopy for thin samples, where the latter even estab- 
lished the existence of individual skyrmions [7]. So far SLPs have been identified in all 
non-centrosymmetric B20 transition metal compounds that order helimagnetically at zero 
magnetic field regardless whether they are pure metals, strongly doped semiconductors or 
even insulators [5, 10, 11]. The excellent theoretical understanding implies that skyrmions 
are a general phenomenon to be expected in a wide range of bulk materials as well as 
nano-scale systems [5, 7, 12], with the identification of spontaneous skyrmion lattices in 
monatomic layers of Fe on an Ir substrate as a major new development [13]. 

Exploratory SANS studies [14] revealed a rotation of the diffraction pattern in the SLP 
of single-crystal MnSi when an electric current applied transverse to the skyrmions exceeded 
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an ultralow threshold of j c ~ 10 6 Am -2 . It is thereby important to emphasize, that j c is 
10 5 times smaller than the currents needed to induce a motion in present-day spintorque 
experiments on ferromagnetic domain walls [15-17]. However, the rotation occurred only in 
the presence of a small temperature gradient (~ 1 K cm -1 ) inducing gradients in the relevant 
forces, which in turn caused rotational torques. It was argued that above the critical current 
density, j c , the skyrmions start to move and that only the sliding skyrmions can be rotated 
by the tiny torques. While microscopic probes such as neutron scattering and Lorentz force 
microscopy may in principle confirm the depinning and motion of the skyrmions, they are 
not capable of detecting the emergent electrodynamics. Instead, to detect the motion of the 
skyrmions, measurements of the emergent electric fields are ideally suited, because they are 
directly proportional to their velocity. This may be readily achieved by means of the Hall 
effect. 

The magnetic properties of MnSi are governed by a combination of ferromagnetic ex- 
change interactions and weak spin-orbit coupling in the absence of inversion symmetry. At 
zero magnetic field MnSi displays a paramagnetic to helimagnetic transition at T c ps 28.5 K. 
In a small range of fields and temperatures below T c the SLP is stabilized [5], where the 
skyrmions (the magnetic whirls) form a hexagonal lattice perpendicular to the applied mag- 
netic field. The lattice constant 2Ah e ii x / v^3 of the skyrmion lattice, determined from the 
reciprocal lattice vectors, is set by the wavelength of the helimagnetic state, A he i ix ~ 180 A. 

Our study of the influence of an electric current on the SLP was carried out on high-purity 
single crystals using a standard six-terminal lock-in technique (see Materials and Methods 
for details). Shown in Fig. 1 (a) are typical temperature dependences of the Hall resistivity, 
p xy , for small currents (black curves). Its dominant features arise from the temperature 
dependence of the anomalous Hall effect. The small maximum in p xy is thereby a charac- 
teristic of the magnetization at the lower boundary of the SLP above a narrow regime of 
phase coexistence between the conical phase and the SLP [18]. Also shown in Fig. 1 (a) is 
the Hall resistivity under an applied DC current density of j = 2.81 • 10 6 Am~ 2 , where a 
suppression (marked in light blue shading) is observed in the temperature range marked by 
the black arrows. The latter is in excellent agreement with the phase boundaries of the SLP 
as inferred from the magnetization and susceptibility. The size of the suppression of the Hall 
signal is similar to the topological Hall contribution Ap xy ps 4nf2cm previously inferred [6] 
from Hall effect measurements at small currents [19]. 
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Detailed Hall data for an applied magnetic field of 0.25 T and a wide range of applied 
DC currents, j, are shown in Fig. 1 (b). From these the evolution of p xy with applied electric 
current density was inferred for selected temperatures as shown in Fig. 2. At temperatures 
outside the SLP, shown in panels (a) and (e), the Hall signal is unchanged as a function 
of j. Within the SLP p xy is within the experimental accuracy unchanged for small current 
densities, j < j c , followed by a clear decrease above j c over a finite range of applied cur- 
rents (light blue shading), which saturates for even larger currents (the method how j c was 
determined is described in the supplement [19]). 

Shown in Fig. 3 (a) is the critical current density as a function of temperature, where 
j c ~ 10 6 A m -2 agrees within a factor of two with the onset of the rotation of the scattering 
pattern observed in the SANS study [14]. As discussed in Ref. [14] the exceptionally low 
value of j c arises from a combination of several factors. First, the very efficient Berry- 
phase coupling between conduction electrons and the spin structure. Second, because of the 
smooth variation of the magnetisation the skyrmion lattice couples only weakly to defects 
and the atomic crystal structure. Third, the long range stiffness and crystalline character 
of the skyrmion lattice [9] leads to a partial cancellation of the pinning forces [21, 24]. 
j c grows by roughly a factor of two when approaching the (first order) transition to the 
weakly field-polarized paramagnetic state. At the same time the suppression of p xy given by 
Ap^ = Ap xy (j ^> j c ) — Ap xy (j <C j c ) is largest in the center of the skyrmion phase with a 
gradual decrease at the lower boundary (cf. Fig. 3 (b)). 

In order to address our observations on p xy from a theoretical point of view we note that 
the forces driving the skyrmion lattice, which also cause the topological contribution to the 
Hall signal, originate in quantum mechanical phases (Berry phases) picked up by electrons 
when their spin follows the orientation n(r,t) = M/\M\ of the local magnetization M. 
The Berry phase can be rewritten [1, 8, 20] as an effective Aharonov-Bohm phase associated 
with 'emergent' magnetic and electric fields, B e and E e . As the Berry phase is given by 
the solid angle covered by h, the emergent fields B e and E e measure the solid angle for an 
infinitesimal loop in space and space-time, respectively 

B i = ^ e ijkh ■ (djh x d k h), E\ = Uh- (<%n x d t n) (1) 

with di = d/dri. Because the sign of the Berry phase depends on the spin orientation, a 
majority spin with magnetization parallel to n carries the emergent charge g| = —1/2 while a 
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minority spin carries the emergent charge g| = +1/2. For a skyrmion, defined as a magnetic 
whirl where h winds once around the unit sphere in the plane perpendicular to B (while n is 
constant in B direction), the total 'emergent flux' is given by f B e dcr = —Arch. It is hence 
topologically quantized to one flux quantum — 2nh/\q e \ per skyrmion (the sign accounts for 
the formation of antiskyrmions in MnSi [5]). Further, according to Eq. (1) and in complete 
analogy to Faraday's law of induction a skyrmion lattice drifting with the velocity v d must 
induce an electric field E e = —v d x B e , where E e /v d inherits its quantization from B e . This 
topological quantization of B e and E e makes skyrmion lattices in metals an ideal system 
to study quantitatively the emergent electrodynamics underlying the coupling of charge and 
magnetism [8]. In MnSi the emergent magnetic field acquires an average strength of 2.5 T, 
i.e. B e ps 2.5T|e/g e |, where e is the electron charge. 

The total force on an electron with momentum k and spin orientation o is therefore given 
[8] by 

F ak = eE + Fh + q%(v akn - v d ) x B e (2) 

where v akn is the velocity of quasi particles in band n, E is the physical electrical field and 
F H <C eE is the Hall force from the normal and anomalous Hall effect. Further dissipative 
drag forces F diss arising for v d ^ are probably much smaller [19] . The extra electric current 
induced by —q%v d x B e transverse to the electrical current has to be canceled exactly by 
the change of the electric Hall field AE± = Ap yx j = —Ap xy j with Ap xy = p xy (j) — Pxy(fy- 
For a current in x and a magnetic field in z direction we find 



Aa yx E Ajj 



E l 
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(v d x B% (3) 

(4) 



where the dimensionless spin polarization P can be obtained by calculating the cross cor- 
relation of the charge current j and the emergent current j e using Kubo formulas. P is 
the ratio of electric currents obtained from E e and E, where a simple approximation for P 
is given in Eq. (4) in the relaxation time approximation of a multiband system (/° CT is the 
Fermi distribution for band n with scattering rate l/r an and spin-orientation a relative to 
the local magnetization). Up to the factor P, the measurement of the Hall field is therefore 
a direct measurement of the emergent field E e and of v d = E e /B e since B e is quantized. 
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The drift velocity v d in absolute units is obtained from Eq. (3) 

e 



v d \\ 



^pin ~ J 



jAp xy _ jA Pxy _ j A fa mm 
B'P pin j c Ap- ^106Am-2 Ap-/- i2 s ' W 
ApZ^ y j c mm 

* 106Am-2°- 12 ~- (6) 



Here we used Ap^ —3 ■ 10~ n f2m in the center of the skyrmion phase and esti- 
mated the effective polarization to be P p=* 0.1 (cf. Ref. [6]). We further used that 
Ap < £ y \e\/(B e P\q e \) pa —v s /j = —\q e /(eM)\P is approximately independent of the local mag- 
netization and therefore of the temperature. As expected, the corresponding pinning ve- 
locities (see labels on the right axes in Fig. 3a) are of the same order of magnitude as the 
electronic drift velocities, thrift ~ j/en ~ 0.16 mm s" 1 for j ~ j c ~ 10 6 Am~ 2 , where we 
estimate n ~ 3.8 ■ 10 22 cm~ 3 from the normal Hall constant [6]. The expression AnhMv v \ n 
may finally be interpreted as the force per skyrmion and per length needed to depin the 
SLP [19]. 

The critical current and the pinning velocities grow by a factor of two just below the 
first order transition T° see Fig. 3(a). This cannot be explained by the decrease of the 
local magnetization (proportional to the small drop in |Ap^|). Instead the increase of j c 
reflects most likely the reduction of the stiffness of the skyrmion lattice when approaching 
Tp, which is only a very weak first order phase transition. With decreasing stiffness the 
skyrmion lattice adjusts better to the disorder potential implying that the pinning forces 
increase substantially [21-23]. Consistent with this picture, there is essentially no tempera- 
ture dependence of j c at the low temperature side of the SLP, where the phase transition is 
strongly first order. Here Ap^ shows a gradual temperature dependence in the range from 
25.8 K to 26.5 K that may be attributed to the phase coexistence between conical phase and 
skyrmion phase as inferred from high precision magnetometry [18]. 

Combining the results of our analysis, we show in Fig. 3 (c) a scaling plot of AE± = 
—Ap xy j, or equivalently of the parallel drift velocity v d \\ and the associated emergent electric 
field E\ as a function of j/j c - To obtain quantitative values for Vd\\ and E\ in physical units, 
one may refer to Fig. 3 (a), Eq. (6), and E\ = Vd\\B e with B e = 2.5T\e/q e \. 

For j < j c the drift velocity vanishes within our experimental precision as the skyrmion 
lattice is pinned by disorder. For j > j c the Magnus forces are sufficiently strong to overcome 
the pinning forces and the skyrmions start to move. For j ^> j c their velocity becomes 
proportional to the current as pinning forces can be neglected. This is precisely the picture 
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which has been developed for the depinning transition of charge density waves and vortices 
[21-24]. However, the precise behavior of skyrmions differs from that of vortices because 
their dynamics is very different. Due to their strong friction, superconducting vortices flow 
approximately perpendicular to the current following the Magnus force and the resulting 
Hall signals are tiny [22, 25]. In contrast, skyrmions drift dominantly parallel to j thus 
reducing the relative speed of spin current and giving rise to a large Faraday field in the 
perpendicular direction. A more technical discussion of the validity of the scaling, the forces 
on the SLP and its direction of motion is given in [19]. 

The direct observation of the emergent electric field of skyrmions reported in this paper 
allowed us to measure their depinning transition and subsequent motion quantitatively. This 
opens the possibility to address fundamental questions of the coupling of magnetism, electric 
currents and defects, respectively. The control and detection of the motion of magnetic 
whirls (skyrmions) by the interplay of emergent and real electrodynamics thereby promises 
to become an important route towards spintronic applications. 
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Methods 



Sample preparation: Single crystals of MnSi were grown by optical float-zoning under 
ultra-high vacuum compatible conditions [26]. The specific heat, susceptibility, and resis- 
tivity of small pieces taken from this single crystal were in excellent agreement with the 
literature, where the residual resistivity ratio was of the order 100. The latter indicates 
good, though not excellent sample purity. The sample quality was the same as for the sam- 
ples studied in the SANS experiments reported in Ref. [14]. Samples for the measurements 
reported here were oriented by Laue x-ray diffraction, cut with a wire saw, and carefully 
polished to size. Current leads were soldered to the small faces of the sample, while Au 
wires for the voltage pick-up were spot-welded onto the surface of the sample. 

Spin torque transport: For our measurements of the Hall and longitudinal resistivity 
we modified a standard six-terminal phase-sensitive detection system such that large DC 
currents could be superimposed on a small AC excitation. The set-up is based on a method 
used for measurements of superconducting tunnel junctions [27]. It was tested on high- 
purity Cu to ensure proper operation. In all experiments the AC excitation amplitude 
were not larger than a few % of the applied DC currents. The samples were carefully 
anchored to the cryogenic system to minimize ohmic heating and temperature gradients. In 
particular, compared to the SANS experiments reported in Ref. [14] all thermal gradients 
were minimized. The Hall signal and the longitudinal resistivity, p xy and p xx , respectively, 
were measured simultaneously at a low excitation frequency of 22.5 Hz. To correct for 
the remaining tiny amount of uniform ohmic heating, which generated a small systematic 
temperature difference between sample and thermometer of less than a few tenths of a K 
for the largest currents applied, we calculated from the longitudinal resistivity p xx the true 
sample temperature. However, we have tested carefully that our results do not depend on 
the precise way how these small ohmic heating effects are corrected. 
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FIG. 1: Temperature dependence of the Hall resistivity in the skyrmion lattice phase of MnSi 
under a large applied DC electric current. To study the effect of the applied DC current it is 
superimposed on a small AC excitation that allows to detect the signal, (a) Hall resistivity for 
various magnetic fields. Under an applied DC current of 2.81 -10 6 Am -2 the Hall signal is suppressed 
in the entire skyrmion phase (blue shading) . (b) Hall resistivity in the SLP for an applied magnetic 
field B = 250 mT under selected applied DC currents between j = and j = 2.81 • 10 6 Am~ 2 . 
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FIG. 2: Typical variation of the Hall resistivity, p X y, of MnSi function of applied DC current 
at B = 250 mT and for selected temperatures. For temperatures in the skyrmion lattice phase, 
shown in panels (b), (c) and (d), the signal is suppressed above an ultralow current density, j c , 
and limits towards a constant value for large currents. Absolute data values are given on the left 
hand side, while the relative change with respect to j = is given on the right hand side. 
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FIG. 3: Main characteristics of the response of the skyrmion lattice phase in MnSi under applied 
electric currents, (a) Temperature dependence of the critical current density, j c . The label on 
the right axis gives an estimate of the pinning velocity, see Eq. (6). (b) Total change of the Hall 
resistivity for large currents, = p xy (j -C j c ) — Pxy(j S> j c ), as a function of temperature, (c) 

Scaling plot of the transverse electric field, AE± = —jAp xy , in units of j c Ap^ induced by the 
moving skyrmion lattice. As AE± is proportional to the emergent field E e , this also constitutes a 
scaling plot of the emergent electric field E e (in units of v p [ n B e ) or of the drift velocity vm in units 
of the pinning velocity t; p i n (shown in (a)). Unsealed data used to construct panel (c) are shown 
in the supplement [19]. 
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Abstract 

Supplementary information for our manuscript, entitled 'Emergent electrodynamics of skyrmions 
in a chiral magnet' is presented. The data are shown in a form essentially as measured, highlighting 
the determination of the critical current density. A detailed technical discussion is presented of 
the forces acting on the skyrmion lattice under electric current flow followed by a discussion of the 
relative size of the Hall effect and the Faraday effect in multi-band and Galilean invariant systems. 



In the following we present at first unsealed data of the main effect reported in the main 
paper. This is followed by a discussion of the forces on the skyrmion lattice. We then 
show how the pinning force and the direction of motion of the skyrmions can be described 
theoretically and present arguments that the dissipative and drag forces on the electrons 
are very small. Finally, we will discuss the relative size of the emergent Hall effect and the 
Faraday effect both in a material with many bands and in a Galilean invariant system. 

I. DRIFT VELOCITIES INFERRED FROM EXPERIMENTAL DATA 

Shown in Fig. 1 are typical data of the product of the current density j with the change 
of Hall resistivity, Ap xy , recorded at an applied field of 250 mT. These form the basis of 
the scaling shown in Fig. 3 (c) of the main text. The data as plotted here was also used to 
determine the critical current density j c at the point where j Ap xy begins to deviate from 
zero. Shown on the right hand side of each panel is the corresponding drift velocity according 
to Eq. (5) of the main text. 

II. FORCES ON THE SKYRMION LATTICE AND DEPINNING TRANSITION 

To understand the origin of the drift velocity v d of the skyrmions, we study the relevant 
forces by projecting the Landau-Lifshitz-Gilbert equation [1] on the translational modes 
using an approach pioneered by Thiele [2]. Besides a Magnus force [3], which is the coun- 
terforce to the emergent Lorentz force arising from B e , and dissipative forces one has also 
to consider pinning forces f p - m arising from inhomogeneities which prohibit the motion of 
skyrmions for small currents. Balancing all forces (per skyrmion, per length, and per KM) 
one obtains 

G x (v s - v d ) + V(f3v s - av d ) = / pin (1) 

/pin ~ -^V pin f(v d /V pin )v d (2) 

where a and (3 parametrize the damping [1], the gyromagnetic coupling vector G = J B e /h 
integrated over the 2d magnetic unit cell evaluates to — AttB, v s = (j e )/M is the velocity of 
the current associated to the charges q e a1 and V = f uc (dh) 2 /2 [2]. Up to tiny corrections of 
the order ~ (fxy/cxx, 

v s is parallel to the applied current. For the skyrmions in MnSi, the 
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FIG. 1: Typical data of the product of the current density j with the change of the Hall resistivity 
Ap xy as a function of current density. The departure of the data points from zero marks the 
critical current density, where the skyrmion lattice begins to move. The scale on the right hand 
side displays the corresponding drift according to Eq. (5) of the main text. 

dissipative terms are expected to be small as a, (5 <C 1 arise mainly from spin-orbit effects 
and spin flip scattering (see Sec. Ill below) and we obtain values of order 1 for T>j\G\ for 
typical parameters [4] where the skyrmion crystal is stable. In the following, we therefore 
neglect the dissipative forces in the parameter range of our experiment. 

The small size of the dissipative forces has important consequences for the motion of the 
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skyrmions, especially in the regime j ^> j c where the pinning force f pin (discussed below) can 
be neglected. In this regime, one obtains v d v s with corrections of order otD/\G\, i.e., the 
skyrmions move parallel to the current in such a way that the Magnus force is canceled. In 
contrast, for vortices in superconductors, where aD/\G\ ^> 1 [5], dissipation dominates such 
that Vd <C v s . Therefore the motion of superconducting vortices is dominantly perpendicular 
to the current following the direction of the Magnus force. Similar to the skyrmions, charge 
density waves [6] flow parallel to the current. However, this is less suprising since the 
dominant force on CDWs is not a Magnus force but arises directly from the electric field. 

The pinning force, / pin , for the moving skyrmion lattice is oriented antiparallel to the drift 
velocity v d - We parametrize its strength by the 'pinning velocity' v v - in and a dimensionless 
function f(vd/v p i n ) with /(0) = 1 and f(x ^> 1) ~ x a , a < 1 (assuming the absence of 
creep). / is not known, but we expect a behavior similar to that of superconducting vortices 
though not identical due to the different dynamics [7-9]. Furthermore, we do not expect 
that / is a fully universal function; it will depend, e.g., on the various elastic constants of 
the skyrmion lattice. Nevertheless, the observed approximate scaling of v d /v pin as a function 
of v s /v pm ~ j/j c shown in Fig. 3 (c) of the main article suggests that / depends only very 
little on the distance to the transition temperature while t> p i n varies strongly. 

Neglecting damping, one can calculate v d /v pin from /. For unkown /, one may determine 
/ from the measurement of v d \\/v pin as a function of j/j c . By projecting Eq. (1) (for a, j3 — 0) 
to the direction perpendicular and parallel to v d and using v\ = v^ + v 2 dl _ and Eq. (5) of 
the main article, one obtains for the components of v d parallel and perpendicular to v s the 
following equations 
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v d \\ 



pm 
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xy 



xy 
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While Eq. (3) (which just copies Eq. (5) of the main article) reflects that the Hall effect 
provides a direct measurement of the drift velocity of skyrmions in parallel direction, the 
other equations show, that - under the assumption that Eq. (1) is valid and that damping 
is negligible - one can also obtain the drift velocity in perpendicular direction, the angle 9 
between and the current and the pinning force as function of Vd/v pin experimentally. 

Unfortunately, our present data is too noisy to extract reliably the force / or v^±. As 
an example, we show in Fig. 2 the dimensionless force / as a function of Vd/v p - m obtained 
from Eqs. (4) and (7). For j < j c the force grows with the applied current but Vd remains 
approximately which results in a vertical line for / < 1, v& — 0. For Vd > 0, the force 
appears to remain approximately constant as a function of Vd, but due to the large noise a 
definite conclusion is presently not possible. 
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FIG. 2: Force / on the skyrmions in dimensionless units as a function of Vd/v p i n from Eqs. (4) 
and (7). 



III. DISSIPATIVE DRAG FORCES AND CORRECTIONS TO THE ADIABATIC 
APPROXIMATIONS 

In this section we present arguments why dissipative drag forces acting on the electrons 
may be neglected. Since we have argued that dissipative forces acting on the skyrmion 
lattice are small this is not surprising. In addition we briefly discuss in this section the role 
of other corrections which may modify the emergent magnetic and electric fields. 
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When the skyrmion lattice starts to move, the electrons experience the forces induced by 
the emergent magnetic and electric fields B e and E e , discussed in the main text. In addition 
dissipative forces F d drag the electrons parallel to the motion of the skyrmion lattice as the 
skyrmion lattice induces effectively also a weak periodic potential arising from at least three 
different sources, notably spin-orbit coupling in the band structure, small variations of the 
amplitude of the magnetization, and variations of B e around its average value. 

To obtain an estimate of the strength of these forces for small drift velocities v d of the 
skyrmion lattice, it is useful to make use of results for the generic problem of a weak time- 
dependent periodic electric potential $(r — v d t) moving in a diffusive system described by 
the diffusion equation for the charge density n 

d t n = V(D(n)Vn) + V(a(n)V$(r - v d t)) (8) 

In a comoving frame of reference one obtains a stationary solution by the transformation 
n — > n(r — v d t) and r — > r + v d t with 

-v d Vn = VD(n)Vn + Va(n)V$(r) (9) 

with the density-dependent diffusion constant D{n) and conductivity a(n). Integrating this 
equation with an integration constant —v d rio and multiplying with v d gives 

- v 2 d (n - n ) = D(n)v d Vn + a(n)v d V$(r) (10) 

Linearizing the equation using the approximation D{n) = D = const, and a{n) = a = const, 
one finds after Fourier transformation with n(r) = J2 q n q e tqr 

nQ * ~ v> + i(Q iVd )D 

where Qi ^ are the reciprocal lattice vectors of the periodic potential. To this linear order 
no dc current is obtained, while only a small periodic modulation of the density is induced 
by the potential. To obtain a uniform electric current, one has to expand one order higher 
in <3>. Expanding D and a to linear order in the changes of n gives rise to an average current 
in the direction of v d 

j = -l/(Vv d ) J(D(n)v d Vn + a(n)v d V<S>(r)) (12) 



dcr^ (v d Qi) 2 $ Qi 2 da ^ $ Q . r 

dn i v* + ((QiV d )D) 2 dn V Q^ D 
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for small v^. For a potential characterized by variations of order AV ~ I^qJ and with 
a ~ ner/m, D ~ v 2 f t, scattering rate 1/r, mean free path /, quasiparticle mass m, Fermi 
energy 6f, Fermi wave length Xp, Fermi momentum hp, and Fermi velocity vf, we can 
therefore estimate the current and the associated dissipative force from 



ej mv d f AV 



2 



(14) 

a t \ e F J 

Note that the current is independent of the scattering rate within our approximation, while 
the dissipative force is proportional to the scattering rate. 

As discussed above, several mechanisms contribute to AV but for a crude estimate of 
its size it is sufficient to note that the change of electronic energy due to AV, which is of 
order (AV) 2 /ep (or much larger), is maximally of the same order as the energy difference 
of the ferromagnetic and the skyrmion state, AE S ~ E M /(kFD) 2 , where E M is a magnetic 
energy scale and D is the radius of skyrmions which for MnSi also describes the distance of 
skyrmions [10]. Therefore we estimate 

/-> ^ hi E M 1 E M 

F d^v d — — v d B e — 15 

D l kfI e_p KfI €f 

For the last estimate, we used that the emergent magnetic field (quantized to one flux 
quantum per magnetic unit cell) is of order of h/D 2 . 

Compared to the force from the emergent electric field, E e = —v d x B e , F d is strongly 
suppressed by two small factors, l/(k F l) and Therefore we expect that the contribution 
of dissipative forces to the topological Hall signal are completely negligible compared to 
the contributions from E e . Only in a tiny regime close to the depinning transition, well 
below our experimental resolution, these contributions might have a chance of becoming 
relevant. In this regime v d is approximately perpendicular to the current and therefore the 
contribution of E e to the Hall effect is suppressed. 

Finally, further corrections to the simple formulas for the effective electric and magnetic 
fields, Eq.(l) in the main text, have to be considered. The formulas for the emergent 
electric and magnetic fields have been derived under the assumption that the conduction 
electron spins follow the local magnetization adiabatically and that one does not have to 
take into account a modification of the bandstructure due to the skyrmion lattice. Due to 
the smoothness of the skyrmion structure and the large distance of skyrmions of the order of 
200A, violation of adiabaticity and such bandstructure effects are probably small suggesting 



that the most important corrections may arise from spin-flip scattering processes which 
scatter electrons, e.g., from a majority to a minority band (for a quantitative discussion 
of corrections to adiabatic transport see, e.g., Ref. [11]). These spin-flip processes, which 
can be interaction- or disorder-induced, arise only due to weak spin-orbit scattering. While 
the non-spin-flip scattering length is estimated to be between 10 and 100 A, the spin-flip 
scattering length is much larger and therefore probably also much larger than the distance 
between the skyrmions ~ 200 A. 

Relative size of emergent Hall- and Faraday effect 

Within the experimental errors, we have found that for rapidly moving skyrmions the 
contribution to the Hall effect arising from the emergent electric field almost cancels the 
contribution from the emergent magnetic field. In this section we argue, that, generically, 
the two different contributions are not identical, while there is an exact cancellation in a 
Galilean invariant system. 

In a Galilean invariant system, a state with finite current is obtained by a Galilei trans- 
formation to a moving coordinate system where skyrmions and electrons drift with the same 
velocity and where therefore the force arising from the emergent magnetic field is exactly 
canceled by the force —q e a [vd x B e ) from the emergent electric field. 

Such a cancellation will not occur in a non-Galilean invariant system and several modifica- 
tions have to be considered. First, the damping constants a and ft in Eq. (1) are generically 
finite and not equal such that the spin velocity v s , which describes the drift of the emer- 
gent current j e , differs from the drift velocity of the skyrmions, Vd, even in the absence of 
pinning. We have, however, argued that damping effects are relatively small, therefore this, 
very likely, will be a small effect and we expect v s v d in a regime where pinning can be 
neglected. As a consequence the average force, {ql(v akn — v s ) x B e ), cancels exactly, but in 
a complex multiband system (with quasiparticle velocity v a k n ) it depends on the orientation 
a of the spin relative to the local magnetization, the momentum k, and the bandindex n) 
one nevertheless expects a finite contribution to the Hall effect. 

For example, this may already be seen in a simple relaxation time approximation to the 
Boltzmann equation assuming for simplicity fc-independent scattering rates, l/r an , which 



8 



depend, however, on the local spin orientation and the band index n. For simplicity, we 
assume a cubic bandstructure without spin-orbit coupling. Within this approximation one 
obtains for o xy two contributions arising from B e and from the emergent electric field for 
large currents where ps v s , describing the emergent Hall and the emergent Faraday effect 

<»*-Xf*^(- BJ t)£i («•) 



m akn \ de J (2tt 

Ef r an q% r-o d 3 k 



2^ y ^.S^^fcnJ { de J (2tt) 3 



where the last fraction describes v s /E = j e /(EM). While both expressions contain certain 
averages over scattering rates, velocities, and the inverse mass l/m fffcn , in a complex multi- 
band system there is no simple relation between cr^J and cr^J. Localized spins, for example, 
contribute only to the term in the denominator of cr^J, thereby reducing the Faraday field. 
The special case of localized spins and one fully spin-polarized band was recently discussed 
in Ref. [12]. But band structure effects can also lead to an enhancement of the emergent 
Faraday effect relative to the Hall effect when M, averaged over all bands, is small. Only 
in a Galilei invarant single-band model with a constant mass the two contributions cancel 
exactly, af y + ofj = 0. 

It is finally interesting to note that under pressure an apparent increase of the topological 
contribution of B e to p xy in MnSi by one order of magnitude has been reported [13]. It is 
presently unclear whether this can be explained by Eq. (16) from a change of the scattering 
rates and/or the band structure with pressure. 
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